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Direct measurement of the cross-shock electric potential at low 
plasma /3, quasi-perpendicular bow shocks 

S. D. Bale\ F. S. Mozer\ and V. V. Krasnoselskikh^ 



We use the Cluster EFW experiment to measure the 
cross-shock electric field at ten low (5, quasi-perpendicular 
supercritical bow shock crossings on March 31, 2001. The 
electric field data are Lorentz-tranformed to a Normal Inci- 
dence frame (NIF) , in which the incoming solar wind veloc- 
ity is aligned with the shock normal. In a boundary normal 
coordinate system, the cross-shock (normal) electric field is 
integrated to obtain the cross shock potential. Using this 
technique, we measure the cross-shock potential at each of 
the four Cluster satellites and using an electric field profile 
averaged between the four satellites. Typical values are in 
the range 500-2500 volts. The cross-shock potential mea- 
surements are compared with the ion kinetic energy change 
across the shock. The cross-shock potential is measured to 
be from 23 to 236% of the ion energy change, with large 
variations between the four Cluster spacecraft at the same 
shock. These results indicate that solar wind fiow through 
the shock is likely to be variable in time and space and re- 
sulting structure of the shock is therefore nonstationary. 



1. Introduction 

At collisionless shocks, a large 'cross-shock' electric po- 
tential arises to oppose the incoming plasma flow. This po- 
tential, and its corresponding electric field, have the sense to 
repel incoming ions (which comprise the bulk of the kinetic 

energy and momentum) and to reflect some ions, providing 
additional dissipation and also plays a role in the redistri- 
bution of the upstream flow and magnetic field to the down- 
stream state; the physics of this energy partitioning is not 
yet fully understood. 

In the Ohm's law sense, the cross-shock electric field arises 
from a combination of the Hall current, electron pressure 
gradients, and drag due to the small population of gyrating 
ions at the shock front. There are few, reliable published 
measurements of the cross-shock potential, largely because 
it requires a DC (double-probe) electric field instrument and 
these have been deployed primarily on low-apogee magneto- 
spheric missions. Formisano [1982] estimated a voltage drop 
of 140 and 240V at two shocks. The double-probe instru- 
ment on ISEE was used to measure electric fields of up to 
100 mV/m at the shock [Wygant et al., 1987] and a Normal 
Incidence Frame (NIF) electric potential of 420V [Scudder 
et al., 1986]. Eastwood et al. [2007] used Cluster measure- 
ments to calculate a potential drop of 260V, which corre- 
sponded to the ion kinetic energy change across the shock. 
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Very large electric fields (600 mV/m), including parallel elec- 
tric fields of 100 mV/m have been measured by the Polar 

spacecraft at a high Mach shock [Bale and Mozer, 2007], 
however the single-spacecraft Polar mission does not allow 
for good estimates of shock frame transformations. Fur- 
thermore, the bow shock is known to have large amplitude 
electric field structure from electron inertial scales [Walker 
et al., 2004] down to Debye scales [Bale et al, 1998]. Sev- 
eral authors have used a Liouville mapping of thermal elec- 
trons across the bow shock to estimate the deHoffman-Teller 
frame shock potential [viz. Schwartz et al., 1988], which can 
be related back to the NIF potential by a frame transforma- 
tion. Gedalin [1997] calculated the expected shock potential 
(assuming pressure balance) and found that it should peak 
at an Alfven Mach number of around 2; this is where the 
magnetic compression begins to saturate, while the shock 
thickness continues to grow [Bale et al., 2003] with Mach 
numbers, leading to a small potential and more magnetic 
reflection at high Mach numbers. 

The four Cluster spacecraft can be used to calculate a 
shock reference frame, which is important to transform cor- 
rectly electric fleld data. In this letter, we measure directly 
the cross-shock electric fleld and compute the NIF electric 
potential. We flnd that the potential varies signiflcantly 
between the different spacecraft at the same shock, indicat- 
ing rapid temporal and/or spatial variations. This seems to 
be consistent with expectations of shock 'reformation' [e.g. 
Krasnoselskikh et al., 2002; Hellinger et al., 2002 Matsukiyo 
and Scholer, 2006] and indicates that plasma flow through 
the shock and particle heating and energization are likely to 
be bursty. 

2. Cluster data 

The four Cluster spacecraft fly together in a controlled 
tetrahedron orbit with apogee near 19 Re and inter- 
spacecraft separations that vary from a few hundred to 
several thousand kilometers. Each winter, apogee passes 
through the dayside of the magnetosphere and Cluster 
crosses the bow shock (at least twice during each 57 hour 
orbit). On March 31, 2001, Cluster encountered the bow 
shock 11 times as a CME/magnetic cloud passed over the 
earth. The large, steady magnetic field of the cloud gives a 
low upstream plasma /3 and hence, very planar bow shocks, 
which are ideal for this study. The alpha particle density was 
especially high during this interval, with an average value of 
rial rip ~ 9%, as is often the case within magnetic clouds. 
The alpha density was used in computing ion masses below; 
however, it is not clear how the enhanced alpha density will 
effect the measured electric potentials. 

The EFW experiment [Gustafsson et al., 1997] measures 
probe-to-spacecraft voltage on four 8 cm spherical voltage 
probes, each extended on wire booms 44 m from the space- 
craft body, in the spin plane. Since only two components 
of the electric field are measured, we assume that E ■ B = 
(ideal MHD) in order to determine the three-component 
electric field and correct this where required later. Since 
the Bz (GSE) component is large here (in a CME), this 
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makes for a good reconstruction of the missing (E^) com- 
ponent of the electric field. The electric field data are cali- 
brated locally near each shock by forcing agreement between 
the component of the ion velocity perpendicular to B and 
{ExB)/B^; this correction is of order 1 mV/mintheX GSE 
direction and minimizes any offsets due to varying plasma 
or photoelectron/secondary electron conditions. The sum 
of the four probe voltages gives an estimate of the space- 
craft floating potential which is related functionally to the 
ambient plasma density. We fit the spacecraft potential lo- 
cally (near each shock) to the plasma density to produce a 
high time resolution 'density proxy' measurement. We also 
use magnetic field data from the PGM experiment [Balogh 
et al., 1997], ion moments from the CIS instrument [Reme 
et al., 1997], and electron temperatures from the PEACE 
instrument [Johnstone et al., 1997]. The fractional alpha 
particle density is estimated using ACE data upstream and 
convected back to the shock crossing time; the alpha den- 
sity is used to estimate the solar wind mass density, where 
needed. 
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Figure 1. The shock coordinate system is defined with 
respect to the normal and the coplanarity plane, which 
contains the magnetic field and velocity vectors. In the 
normal incidence frame, the solar wind velocity is di- 
rected along the normal (red arrow). The convection 
electric field E — —v x B (green arrow) is normal to the 
coplanarity plane and is approximately constant through 
the shock layer. 



3. IMeasurements at bow shock crossings 

Of the 11 bow shock crossings observed by Cluster on 
March 31, 2001, ten (10) of them were chosen for this study. 
These shocks were used by Maksimovic et al. [2003] to study 
global shock motion. 

3.1. Shock normals 

We compute the shock normal h by comparing the shock 
arrival time at the four Cluster satellites Ti and inverting the 
matrix equation r ■ n = Vsht, where r is a matrix of relative 
spacecraft positions and Vsh is the shock speed in the space- 
craft frame [viz Bale et al., 2003]. We do this separately 



using both spacecraft potential and magnetic field magni- 
tude data and find normals that agree to within a few de- 
grees and speeds that agree to within a few km/s, typically. 
Minimum variance normals are also in good agreement. A 
single spacecraft crosses the shock in 2-10 seconds typically 
and the shock transit time between different spacecraft is 
from 7 to 37 seconds. 
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Figure 2. A typical shock crossing on March 31, 2001. 
The electric and magnetic fields and perpendicular ve- 
locity are in the coordinate system described in Pigure 1; 
the purple line is the x component, green is y, and red is 
the z component, large dots in the bottom panel show the 
perpendicular component of the ion velocity measured by 
CIS. The electric potential (top panel) is computed by in- 
tegrating the normal (purple) component of the electric 
field (middle panel). The By (out-of-coplanarity) mag- 
netic field can be seen in the second panel (green trace). 



3.2. Lorentz frame and coordinate transformations 

The spacecraft-frame (measured) electric field can differ 
from that of the shock frame by several mV/m. There- 
fore, we Lorentz transform the measured field into the shock 
frame using the measured shock velocity Vah = Vshn to gen- 
erate the shock-frame electric field Esh = E + Vsh x B. Then 
we compute the NIP velocity ujv/f — hx (vu x n), where Vu 
is the solar wind velocity in the shock frame, and Lorentz 
transform the electric field Enif = Esh + vnif x B and the 
incoming flow velocity Vsw ~ Vu ~ vnif to the NI frame. 
Note that to this point, we have assumed that E ■ B = 
(which is a Lorentz invariant). 

Finally, the shock frame electric and magnetic flelds and 
velocity are rotated into a coordinate system (Pigure 1) 
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which is defined with the shock normal as the x direction 
and the maximum variance of the magnetic field vector as 
the z direction; x and z define the coplanarity plane. In this 
coordinate system, the cross-shock electric field is in the x 
direction, the convection electric field — Wsm x _B is in the 
y direction, and the magnetic field shearing direction is z. 
Figure 2 shows data in this frame and coordinate system at 
one of our shock crossings. 

3.3. The cross-shock electric potential 

The shock speed calculated as described in Section 3.1 is 
used to generate a spatial shock profile dx — Vshdt and the 
cross-shock {x component) electric field can be integrated 
directly to obtain the cross-shock potential. To compensate 
for the assumption of E ■ B = 0, we now assume that the 
true cross-shock electric field lies purely in the x (normal) 
direction, and that we are measuring only the projection of 
it perpendicular to 13; therefore the cross-shock field can be 
written as iSji = E^x = {Es smQtn)-i-~ {Eg cosQtn)\\, where 

II and J- are unit vectors parallel and perpendicular to the 
magnetic field (and hence in the coplanarity plane). The 
perpendicular component (underlined above) is the mea- 
sured E±^x and therefore we can recover the cross-shock field 
amplitude as Es — E±^x/ sin0i,„ and this is the field that 
we integrate to obtain the NIF potential; this represent a 
correction of from 0.1% for 87° to 21% for G^n 56° 

(see Table 1). It is interesting to note that the cross-shock 
electric field is comprised of both a large-scale DC field and 
shorter wavelength, spiky structures of comparable ampli- 
tude (viz [Walker et al, 2004]). These structures are in- 
cluded in the integral of cross-shock potential, however, the 
resulting potential profile is relatively smooth and mono- 
tonic (until the magnetic overshoot, which is mimicked in 
the electric potential profile). 

At each shock, the electric potential 4)sci is computed, 
as described above, on each of the four Cluster spacecraft 
{i = 1, 4) and an 'average' potential (j)A is computed by first 
aligning (in time) the data from the four spacecraft, com- 
puting an average cross-shock electric field profile, and then 
integrating it, so that (pA is like a potential of the volume- 
averaged field. Table 1 lists the 10 shocks, their macro- 
scopic parameters, the Alfven Mach number Ma, electron 
and proton plasma beta, the shock tangent angle Obn, and 
the ion energy E = l/2mv'^ and its change across the shock 
AE = l/2m(w„ — w|), along with the measured electric po- 
tentials. 



4. Discussion 

The last columns of Table 1 show (jtA/E and <j>A/ AE, the 
average cross-shock potential normalized to the upstream 
ion kinetic energy E = lj2mvt and its change across the 
shock AE = l/2m{vl-vi). This is a measure of the ability 
of the shock to oppose the directed plasma fiow; i.e. when 
the electric potential (fiA approaches the upstream ion en- 
ergy E the shock should turn back the entire solar wind 
thermal ion population. Figure 3 shows (j>/ E with the black 
dots as 4)A and the error bars representing the maximum 
and minimum at each shock and plotted against Alfven 
Mach number; red dots show the values of (\>aI AE The elec- 
tric potential can vary by nearly 100% in some cases and in 
three of the ten cases, (j)A/ E is greater than 1. 
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Figure 3. The cross-shock potential normalized to 
the upstream NIF ion kinetic energy E = (t>/{l/2mv1) 
plotted against Alfven Mach number. The black dots 
are computed from the 'average' shock potentials (j)A 
and the error bars show the range of potential varia- 
tion between the 4 Cluster satellites. Red dots show 
the ratio of (j)A to the NIF ion kinetic energy change 
AE = l/2m{v^ — v^). The dotted line is the analyti- 
cal relationship from Gedalm (1997). 



Table 1. Normal Incidence Frame (NIF) shock parameters and measured cross-shock electric potentials on March 31, 2001. 



cross-shock potential (volts) 



Shock Time 


Ma 


/9e 






E (cV)'' 


AE (eV)*^ 




4>SC2 


4>sc-i 


<f>SC4 




<t>AlE 


<i>AlAE 


17:14:45 


2.4 


0.03 


0.02 


83° 


2299 


1926 


1973 


1900 


2057 


3245 


2190 


0.95 


1.14 


17:18:50 


2.9 


0.01 


0.01 


85° 


2223 


2116 


1200 


830 


1754 


1232 


1223 


0.55 


0.58 


17:36:47 


3.2 


0.06 


0.05 


86° 


2458 


2110 


251 


535 


825 


755 


482 


0.20 


0.23 


17:38:20 


3.9 


0.05 


0.06 


86° 


2561 


2331 


840 


1884 


923 


1375 


1239 


0.48 


0.53 


18:02:15 


3.4 


0.03 


0.03 


87° 


2818 


2467 


893 


714 


1126 


971 


909 


0.32 


0.37 


18:28:40 


5.5 


0.10 


0.10 


84° 


2487 


2266 


2373 


1541 


1185 


2348 


1800 


0.72 


0.79 


18:48:20 


2.5 


0.11 


0.07 


57° 


2581 


2121 


1748 


1127 


920 


1146 


1039 


0.40 


0.49 


19:00:41 


3.7 


0.10 


0.09 


64° 


2362 


2053 


2968 


2813 


3623 


2992 


2791 


1.18 


1.36 


19:46:37 


2.6 


0.02 


0.02 


62° 


2261 


1798 


1012 


809 


726 


980 


799 


0.35 


0.44 


21:34:06 


2.7 


0.03 


0.03 


56° 


1820 


1688 


4303 


4987 


3495 


5402 


3992 


2.19 


2.36 



^ E = l/2mv^ is the NIF upstream ion kinetic energy 

^ AE = l/2m{vl - vj) is the NIF ion kinetic energy change from upstream to downstream 
Note that ipA is not the average of the four potentials, rather it is the cross-shock potential of the average shock electric field 
profile. 
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Gedalin [1997] estimated the cross-shock potential in the 
NIF analytically (assuming a monotonic shock profile and 
pressure laalance) and found that the cross-shock potential 
peaks at small Alfven Mach numbers [Ma ~ 2), qualita- 
tively in agreement with dHT potentials inferred from Liou- 
ville mapping of thermal electrons [Schwartz et al.,, 1988]. 
Our NIF potentials show a similar trend (Figure 3), al- 
beit with poor statistics; the dHT and NIF cross-shock 
{x) electric fields are related by a frame transformation 

EdHT = Enif + (Ey / Bx)By. 

Theoretical studies of shock reformation or nonstationar- 
ity predict that quasiperpcndicular shock fronts should be 
unstable in certain regimes of Mach rmniber and plasma (5 
[e.g. Krasnoselskikh et al, 2002; Hellinger et al, 2002; Mat- 
sukiyo and Scholer, 2006] . In particular, the constraint that 
the incoming solar wind speed be larger than the whistler 
phase speed and/or the ion thermal speeds allows the de- 
velopment of a shock front instability that results in 'over- 
turning' of the front with a characteristic timescale of the 
ion gyroperiod and a spatial scale of the ion gyroradius. 

Recently Lobzin et al. [2007], using Cluster observa- 
tions, provided convincing evidence that high Mach num- 
ber quasiperpcndicular shocks arc nonstationary, moreover, 
a quasi-periodic shock front reformation takes place which 
modulates the reflected ion population. 

In the reformation scenario, the shock electric potential 
will oscillate to large values on timescalcs of the ion gyrope- 
riod Tci- The time between shock crossings here (from space- 
craft to spacecraft) ranges from 7 to 30 seconds, while the 
ion cyclotron period is Td ~ 2 s; so it is plausible that these 
shocks are reforming rapidly and the four Cluster space- 
craft each encounter the same shock at a different phase 
of the reformation cycle resulting in large variations in the 
measured potential. This strongly varying electric potential 
should produce a modulated reflected ion flux, as observed 
by Lobzin et al, [2007]. During this interval, the Cluster 
spacecraft were separated by 400-900 km, while pi ~ 200 
km, so that the spacccraft-to-spacecraft variations may be 
spatial, rather than temporal. 

This is the first multi-spacecraft study of the cross-shock 
electric potential (to our knowledge) and it shows that shock 
electric potentials vary largely on the timescale of the ion 
gyroperiod and/or spatial scales of ion gyroradii and that 
the Normal Incidence frame potential is often observed to 
be greater than the ion kinetic energy change across the 
shock. These observations are consistent with the quasiper- 
pcndicular shock reformation scenario and suggest that the 
transmission and reflection at the shock front arc a bursty 
phenomena. This shock reformation may be due to an inher- 
ent instability (as described in references above) or due to 
solar wind driving, although there is no one-to-one signature 
of this in the upstream data [Maksimovic et al., 2003]. 

Acknowledgments. Cluster EFW data analysis at SSL is 
supported by NASA grant NNG05GL27G to the University of 
California. We acknowledge the Cluster CIS, PEACE, and FGM 
and the ACE EPAM teams for data. 



References 



tion region: evidence of electron phase space holes, Geophys. 
Res. Lett., 25, 2929, 1998. 
Bale, S. D., F. S. Mozcr, and T. S. Horbury, Density-transition 
scale at quasi-perpendicular coUisionless shocks, Phys. Rev. 
Lett., 91, 265004, 2003. 

Bale, S. D. and F. S. Mozer, Measurement of large parallel and 
perpendicular electric fields on electron spatial scales in the 
terrestrial bow shock, Phys. Rev. Lett., 98, 20501, 2007. 

Balogh, A. et al.. The Cluster Magnetic Field Investigation, Space 
Science Rev., 79, 65, 1997. 

Eastwood, J. P., S. D. Bale, F. S. Mozer, and A. J. Hull, Contri- 
butions to the cross shock eleetrie field at a quasiperpendicular 
coUisionless shock, Geophys. Res. Lett., 34, L17104, 2007. 

Formisano, V., Measurement of the potential drop across the 
Earth's coUisionless bow shock, Geophys. Res. Lett, 9, 1033, 
1982. 

Gedalin, M., Ion heating in oblique low-Mach number shocks, 

Geophys. Res. Lett, 24 , 2511, 1997. 
Gustafsson, G., et al., The Electric Field and Wave experiment 

for the Cluster mission, Space Science Rev., 79, 137, 1997. 
Hellinger, P., P. Travnicek, and H. Matsumoto, Reformation of 

perpendicular shocks: hybrid simulations, Geophys. Res. Lett, 

29, 2234, 2002. 

Johnstone, A. D., et al., PEACE: A Plasma Electron And Current 
Experiment, Space Science Rev., 79, 1572, 1997. 

Krasnoselkskikh, V. V., B. Lembege, P. Savoini, and V. V. 
Lobzin, Nonstationarity of strong coUisionless quasiperpendic- 
ular shocks: theory and full particle simulations, Phys. Plas- 
mas, 9, 1192, 2002. 

Lobzin, V. V., V. V. Krasnoselskikh, ,I.-M. Bosqued, ,I.-L. Pineon, 
S. .1. Schwartz, and M. Dunlop, Nonstationarity and reforma- 
tion of high-Maeh-number quasiperpendicular shocks: Cluster 
observations, Geophys. Res. Lett., 34, L05107, 2007. 

Maksimovic, M., S. D. Bale, T. S. Horbury, and M. Andre, Bow 
shock motion observed with Cluster, Geophys. Res. Lett, 30, 
1393, 2003. 

Matsukiyo, S. and M. Scholer, On reformation of quasi- 
perpendicular coUisionless shocks, Advances Space Res., 38, 
57, 2006. 

Reme, H. et al., The Cluster Ion Spectrometry (CIS) experiment, 

Space Science Rev., 79, 303, 1997. 

Schwartz, S. J., M. F. Tliomsen, S. J Bame, and J. Stansberry, 
Electron heating and the potential jump across fast mode 
shocks, J. Geophys. Res., 88, 12923, 1988. 

Scudder, J. D., A. Mangeney, C. Lacombe, C. C. Harvey, and 
T. L. Aggson, The resolved layer of a coUisionless, high /3, su- 
percritical, quasi-perpendicular shock wave 2. Dissipative fluid 
electrodynamics, J. Geophys. Res., 91, 11,053, 1986. 

Walker, S. N., H. St. C. K. AUeyne, M. A. BaUkhin, M. Andre, 
and T. S. Horbury, Electric field scales at quasi-perpendicular 
shocks, Annates Geophys., 22, 2291, 2004. 

Wygant, ,1. R., M. Bensadoun, and F. S. Mozer, Electric field 
measurements at suberitieal, oblique bow shock crossings, J. 
Geophys. Res., 92, 11109, 1987. 



Bale, S. D., P. J. Kellogg, D. E. Larson, R. P. Lin, K. Goetz, and 
R. P. Lin, Bipolar electrostatic structures in the shock transi- 



Stuart D. Bale, Space Sciences Laboratory, University of Cal- 
ifornia, Berkeley, CA 94720-7450, USA. (bale@ssl.berkeley.edu) 



